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METHOD AND APPARATUS FOR PARAMETER 
ESTIMATION IN A GENERALIZED RAKE RECEIVER 

BACKGROUND OF THE INVENTION 
[0001] The present invention generally relates to wireless communication services, and 
particularly relates to tracking signal impairment correlations of received communication signals. 
[0002] RAKE receivers are well known in the communication arts and find widespread use 
in Code Division Multiple Access (CDMA) systems, such as in IS-95, IS-2000 (cdma2000), and 
Wideband CDMA (WCDMA) wireless communication networks. The name derives from the 
rake-like appearance of such receivers, wherein multiple, parallel receiver fingers are used to 
receive multiple signal images in a received multipath signal. By coherently combining the 
finger outputs in a RAKE combiner, the conventional RAKE receiver can use multipath reception 
to improve the Signal-to-Noise Ratio (SNR) of the received multipath signal. 
[0003] However, as is known to those skilled in the art, the conventional RAKE receiver is 
optimal only in certain limited circumstances. For example, the presence of self-interference 
and multi-user access interference both degrade the performance of a conventional RAKE 
receiver. To that end, the assignee of the instant application has made application for one or 
more patents relating to the use of a "generalized" RAKE receiver architecture, wherein receiver 
performance is improved by increasing the sophistication of combining weight generation. 
[0004] Thus, in the generalized RAKE architecture, the combining weight calculations 
consider correlations of one or more signal impairments across RAKE fingers. For example, a 
generalized RAKE receiver may track noise correlations across those fingers. Generalized 
RAKE receivers also may include a comparatively larger number of fingers such that extra 
fingers may be positioned off of the signal path delays. Indeed, a generalized RAKE receiver 
can gain performance improvements by shifting these extra fingers to optimize the SNR of the 
received signal. Correlations of signal impairments can also be used in SNR estimating often 
referred to as signal to interference ratio (SIR) estimation. SIR estimation is used in power 
control as well as in monitoring link quality and rate adaptation. 
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[0005] Regardless, by using its knowledge of how selected signal impairments are 

correlated across fingers, the generalized RAKE receiver can compensate the finger combining 

weights such that receiver performance is improved. Of course, the need to determine signal 

impairment correlations with sufficient accuracy and rapidity stands as a primary challenge 

associated with implementation of the generalized RAKE receiver. 

SUMMARY OF THE INVENTION 
[0006] The present invention provides a method and apparatus to estimate signal 
impairment correlations for one or more received signals of interest using a model-based 
technique wherein the model is adapted responsive to recurring measurements of signal 
impairment correlations that can be made on a frequent basis, e.g., every timeslot of a 
Wideband CDMA (WCDMA) frame, thereby dynamically tracking even rapidly changing signal 
impairment correlations. In one or more exemplary embodiments, the present invention thus 
comprises a method of determining received signal impairment correlations for use in generating 
RAKE combining weights and/or SIR estimates, wherein the method comprises providing a 
model of received signal impairment correlations comprising one or more impairment terms 
scaled by corresponding model fitting parameters, and adapting each of the model fitting 
parameters responsive to recurring measurements of the received signal impairment 
correlations such that the model of received signal impairment correlations dynamically tracks 
changing reception conditions. 

[0007] In at least one embodiment, an exemplary method comprises providing an 
interference impairment term scaled by a first fitting parameter and a noise impairment term 
scaled by a second fitting parameter to model received signal impairment correlations, 
measuring received signal impairment correlations at each of one or more successive time 
instants and, at each time instant, fitting the model to measured received signal impairment 
correlations by adapting instantaneous values of the first and second fitting parameters, and 
maintaining the model by updating the first and second fitting parameters based on the 
instantaneous values calculated for them at each time instant. 
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[0008] Embodying one or more of these exemplary methods, an exemplary wireless 
communication terminal for use in a wireless communication network comprises a radio front- 
end circuit, e.g., a radio processor, configured to provide one or more received signals of 
interest corresponding to one or more antenna-received signals, and a RAKE receiver circuit 
configured to generate one or more RAKE combined signals by RAKE processing the one or 
more received signals of interest. An exemplary RAKE receiver circuit is configured to calculate 
RAKE combining weights and SIR estimates by providing a model of received signal impairment 
correlations for a received signal of interest comprising an interference impairment term scaled 
by a first fitting parameter and a noise impairment term scaled by a second fitting parameter, 
measuring received signal impairment correlations at each of one or more successive time 
instants and, at each time instant, fitting the model to measured received signal impairment 
correlations by adapting instantaneous values of the first and second fitting parameters, and 
maintaining the model by updating the first and second fitting parameters based on the 
instantaneous values calculated for them at each time instant. 

[0009] Of course, the present invention includes additional features and advantages as 
highlighted in the following detailed discussion. Those skilled in the art will recognize additional 
features and advantages upon reading that discussion, and upon viewing the accompanying 
drawings. 

BRIEF DESCRIPTION OF THE DRAWINGS 
[0010] Fig. 1 is a diagram of an exemplary method of modeling received signal impairment 
correlations according to the present invention. 

Fig. 2 is a more detailed diagram of an exemplary method of modeling received signal 
impairment correlations according to the present invention. 

Fig. 3 is another exemplary embodiment of modeling received signal impairment 
correlations. 
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Fig. 4 is a diagram of a simplified transmitter and receiver, wherein the receiver includes 

a generalized RAKE receiver circuit according to one or more embodiments of the present 

invention. 

Fig. 5 is a diagram of the generalized RAKE receiver circuit of Fig. 4. 
Fig. 6 is a diagram of an exemplary combining weight and SIR generator included in the 
circuit of Fig. 5. 

Fig. 7 is a diagram of an exemplary wireless communication network, in which a mobile 
terminal includes a receiver circuit configured according to the present invention. 

DETAILED DESCRIPTION OF THE INVENTION 
[0011] In accordance with one or more exemplary embodiments of the present invention, a 
receiver, such as the receiver circuit of a mobile terminal for use in a wireless communication 
network, includes a generalized RAKE receiver circuit that generates RAKE combining weights 
and SIR estimates based on modeled received signal impairment correlations. As is known in 
the RAKE receiver arts, the output signals obtained from the individual RAKE fingers can 
include cross-correlated "impairments." As used herein, the term "impairment" has a broad 
definition that includes, but is not limited to, one or more of the following items: self and multi- 
user interference and noise. For additional background regarding exemplary generalized RAKE 
operations, one may refer to the co-pending and commonly assigned U.S. patent application 
assigned Ser. No. 09/344,899, which is incorporated herein by reference. 
[0012] In the generalized RAKE architecture, denoted as "G-RAKE" herein, the combining 
weights used to combine the RAKE finger outputs are compensated for cross-finger correlations 
in one or more of these impairment terms, and the present invention provides a method and 
apparatus for carrying out such compensation using a model-based approach. Fig. 1 broadly 
outlines an exemplary method to implement the present invention's model-based method of 
impairment correlation estimations. In the context of Fig. 1 , the exemplary model comprises one 
or more impairment terms, which may comprise structured elements, e.g., particularly defined 
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matrices, with each term being scaled by a model fitting parameter. This model structure is 

explored in much greater detail later herein. 

[0013] In examining the illustrated processing logic, one sees that processing "begins" with 
providing a model of received signal impairment correlations (Step 100). Those skilled in the art 
will appreciate that "providing" in this context may comprise configuring a Digital Signal 
Processor (DSP), microprocessor, or other processing logic, to implement an impairment 
correlation model according to stored computer instructions. With the model available, the 
exemplary receiver estimates impairment correlations from received signal measurements at 
successive instants in time, and then updates the model based on these measurements. In this 
sense, then, the present invention provides an adaptive modeling method that "learns" 
impairment correlations directly from received signal measurements. 

[0014] Thus, if it is time to make such measurements and update the model (Step 102), the 
exemplary receiver logic measures received signal impairment correlations — at least for the 
impairment terms of interest — (Step 104), and calculates instantaneous values for the model's 
fitting parameters (Step 106). These instantaneous values are then used to update the model's 
fitting parameters, such that the model of received signal impairment correlations tracks 
changing reception conditions (Step 108). 

[0015] Turning to exemplary model details, the impairment correlation matrix R of a received 
signal for a mobile terminal receiving CDMA signal transmissions from a radio base station can 
be expressed as a function of certain parameters based on the theoretical expression given as, 

R=^-R, 0) 

where E p is the pilot energy per unit time, E t is the total base station energy per unit time, N is 
the spreading factor, C is a scaling factor, N 0 is the noise factor, R/ is an interference correlation 
matrix, and R„ is a thermal noise correlation matrix arising from the autocorrelation properties of 
receiver filtering. Note that R 7 may be constructed as an interference covariance matrix and R„ 
may be constructed as a noise correlation matrix. The terms "covariance" and "correlation" as 
used herein should be understood as interchangeable unless the context of a particular passage 
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makes an explicit distinction between the two terms. Those skilled in the art will, of course, 

appreciate that covariance is a special case of cross-correlation with zero mean. 

[0016] In looking at Eq. (1) above, it should be noted that a receiver generally cannot know 

E t IE p and N 0 explicitly. One or more exemplary embodiments of the present invention 

obviates this problem by determining the model's impairment terms (R 7 and R„ ) explicitly using 

channel coefficient and receiver pulse shape information. Those skilled in the art will appreciate 
that a given receiver can be configured with knowledge of its receiver filter pulse shape, e.g., its 
filter autocorrelation function, and can maintain channel coefficient estimates based on receiving 
pilot symbols, training data, or other signal known a priori to the receiver such that reception of 
the known signal can be used to characterize the propagation channels. Exemplary formulas 
are given herein for computing R 7 and R„ in terms of channel coefficient and pulse shape 
information. Thus, the present invention provides an impairment correlation model that, in an 
exemplary embodiment, comprises an interference term that is scaled by a first model fitting 
parameter, and a noise term that is scaled by a second model fitting parameter. Using this 
method, the impairment correlations, R, can be modeled as, 

R = aR /+ pR„ (2) 

where 

R/k^ 2 )=ZzW/ Y R^-m^-OR/fe-^-rJ 

1=0 q=0 m=~o,ro*0 (3) 

R.(rf 1 ,rf 2 ) = R J ,(rf 1 -rf 2 ) 

and where g{r)=j!g l S(T-T l ) (complex channel model, pilot channel), R 0 (r) is a pulse shape 

autocorrelation function, T c is a CDMA chip period, and d k is the delay of the kth G-RAKE finger. 
Note that the g values are channel coefficients corresponding to the pilot channel, i.e., channel 
coefficients estimated directly from a received pilot channel signal. 

[0017] An exemplary model-based impairment estimation method based on the above 
equations comprises a couple of basic steps. First, the impairment model is adapted based on 
received signal impairment correlations estimated from measurements taken at successive time 
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instants. At each time instant, the method provides an instantaneous estimate of the fitting 

parameters a and p , while a second step may or may not smooth the estimates to eliminate 

estimation noise and provide better receiver performance. 

[0018] An exemplary method according to the above information is performed at each of a 
number of repeating time intervals, e.g., slots, and comprises: 

a. Measure impairment correlations for the received signal in the current slot, i.e., 
determine a rough estimate of impairment correlations expressed as matrix 

R(slot) ; 

b. Calculate per-slot model terms R f (slot) and R n (slot) ; 

c. Determine instantaneous model fitting parameters a inst and J3 inst for the slot 
based on performing a Least Squares fit of, 

R(slot) m a inst Rj (slot) + p inst R n (slot) ; (4) 

d. Update the model fitting parameters a and fj based on the instantaneous fitting 
parameters, e.g., update filtered values of a and/? using a inst and p inst ; and 

e. Calculate the modeled impairment correlations R(slot) to be used in generating 
RAKE combining weights and an SIR estimate for the current slot as, 
R(slot) = cRj(slot) + j3R n (slot) . (5) 

Note that R n (slot) may be denoted simply as R n since the thermal noise characteristics typically 
do not change significantly from slot to slot and thus can be updated more slowly. Also, it 
should be appreciated by those skilled in the art that the term slot is given broad construction 
herein, and can denote any time period, such as a regularly repeating frame time in a wireless 
communication signal. Still further, those skilled in the art will appreciate that parameter 
updating may be performed at irregular intervals and/or as needed, such as in response to the 
SNR after combining falling below a defined threshold. 

[0019] According to the above exemplary method, the first step in the model-based 
impairment correlation estimation process is based on generating a rough — e.g., potentially 
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noisy — estimate of received signal impairment correlations based on measurements of the 
received signal. Thus, the exemplary receiver is configured to generate measurements of the 
impairment correlation matrix on a short-term basis (every slot, every other slot, etc.). Given this 
measurement, knowledge of the channel (coefficients g i and delays r/), and the RAKE finger 
delays d, every quantity in Eq. (4) is known except for a inst and J3 insr One can isolate these 

unknowns by rewriting Eq. (4) in an equivalent form based on stacking of the columns of R , 
R 7 and R„, which yields, 







r /,l r ",l 
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r /,2 r »,2 
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inst 
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inst 



(6) 



where ?/ = I th column of measured impairment correlation matrix R , 17 , = I th column of 

interference correlation matrix R 7 , and r n i - i th column of thermal noise correlation matrix R„. 

(With zero means, the correlation matrices may be understood to be covariance matrices.) 
[0020] As specified above, this system of equations can be solved via least squares. An 
exemplary least squares solution is given as, 



inst 



\_Pinst 



(7) 
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where 



P = 
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[0021] Note that this least square fitting need not be applied to all elements in the measured 

impairment correlation matrix R . For example, only the diagonal elements and the first off- 
diagonal elements could be used to perform the fitting operation. In general, the impairment 
correlation matrix is Hermitian symmetric, so only the diagonal and one of the upper or lower 
triangles would be used. This property can be used throughout, so that only unique matrix 
elements need to be computed and stored by the exemplary receiver. Also, it may be desirable 
to restrict the instantaneous model fitting parameter estimates that result. For example, any 
estimates less than zero could be reset to zero and the receiver could be configured to then 
solve for the other estimate(s). 

[0022] After obtaining the instantaneous fitting parameter values, the exemplary receiver 
uses them to update the long-term model fitting parameters. The exemplary receiver can be 
configured with essentially any low-pass filter to smooth the instantaneous fitting parameter 
estimates. One efficient filter is given by 



where n denotes the current slot, and where 0 < ^ *• 

[0023] Fig. 2 sets out the above method's operations in exemplary processing logic that can 
be implemented in the receiver using an appropriately programmed DSP, microprocessor, or the 
like. The illustrated logic assumes that the receiver is configured to provide a set of combining 
weights for every slot, such as for every traffic channel slot in a WCDMA system, wherein the 
receiver receives ten (10) pilot channel symbols on a Common Pilot Channel (C-PICH) signal 
during each slot. 

[0024] For each slot the receiver obtains desired (despread) pilot (CPICH) symbols and 
RAKE finger delays (Step 110), and estimates net and medium channel coefficients (Step 112) 
using the CPICH, as follows. First, for only fingers that correspond to paths (channel delays), 
we compute 




(10) 
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where h comprises net channel coefficients, x(i) is a vector of despread CPICH symbols 

corresponding to channel delays for symbol index i, is a reference CPICH symbol for 

symbol index i. In general, embedded pilot symbols or data symbols and decision feedback 
may also be used. 

[0025] Next, the exemplary receiver estimates medium coefficients for channel delays 
based on, 



g = 



-1-1 



(12) 



R p {*l-\ -To)R p (t L -\ -*\) - R p (rL-\ -ft-i). 
and then estimates net channel coefficients for all the G-RAKE fingers via 

~R p {d 0 -x 0 ) R p (d 0 -x x ) ••• R p {d 0 -x L _ x )' 
R p {d x -x 0 ) R p {d x -x x ) - R p {d x -x L _ x ) 

R p (d K -x 0 ) R p {d K -x x ) - R p {d K -x L _ x ) 

[0026] The exemplary receiver then computes the measured impairment correlation matrix 
(Step 114) as, 



h = 



(13) 



R=^i(x(«V0)-ii)(xOVG)-h) w , 

y i=0 



(14) 



and computes the model's impairment correlation terms — interference and noise — (Step 116) 
based on, 



L-\L-\ 



^/(^l^2)=ZZ^* Z R p i d \ - mT c - *l K *(di-mT c -T q ) 
l=0q=0 m=oo,m*0 

*n(d x ,d 2 ) = R p {d l -d 2 ) 



(15) 



Note that in practice, the infinite summations may be truncated. The exemplary receiver then 
optionally calculates instantaneous values of the model fitting parameters a ins t and p^f (Step 
118) using Eq. (7) given above. (It should be understood that smoothing may be skipped for 
one or both of the interference impairment matrix scaling factor and the noise impairment matrix 
scaling factor.) 
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[0027] Processing continues with the receiver updating the model fitting parameters a and p 
based on the instantaneous values using Eq. (10) above (Step 120). The receiver then 
computes the modeled impairment correlation matrix, R(slot) , to be used in generating the 
RAKE combining weights, and an SIR estimate for the current slot using the results from the 
above steps (Step 122). The receiver then can compute a SIR estimate and the RAKE 
combining weights used by the G-RAKE receiver circuit such that the combination of the 
different signal images is carried out in consideration of the modeled impairment correlations 
between those images (Step 124). These steps generally are repeated for each of a series of 
successive slots (Step 126). 

[0028] The G-RAKE receiver forms combining weights using R(slot) , and h by basically 
computing: 

w = R- l (slot)h . (16) 
Rather than inverting a matrix, an iterative approach such as Gauss-Seidel can be used. The 
decision variable or symbol estimate is obtained by combining the traffic despread values, 
denoted y, giving 

z = w"y. (17) 

The G-RAKE receiver also forms a SIR estimate using R(slot) and h by basically computing: 

SIR = h^R" 1 {slotjh = h" w (18) 
where superscript H denotes Hermitian transpose. 

[0029] As for initialization of the impairment correlation model, one approach would be to 
initialize using the RAKE solution. This involves setting or to 0 and fi to a positive value, such as 
1 or an estimate of noise power, obtained by conventional means. Of course, other forms of 
initialization are possible and may be used as needed or desired. 

[0030] While specific approaches for channel estimation and parameter fitting are given, the 
invention is not limited to these particular approaches. Below is described a second 
embodiment which employs different approaches for channel estimation and parameter fitting. 
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This embodiment offers an exemplary method of computing channel estimates based on an 
interpolate approach. The summation limits in the below equations assume a WCDMA slot 
wherein one transmit antenna sends ten pilot symbols per traffic slot. 
[0031] With the second embodiment, the medium coefficients are set equal to the net 
coefficients (g = h). The net coefficients are computed for all fingers as an interpolative function 
of time within the slot. Specifically, Step 1 12 computes the net and medium coefficients using 

hW=[^]i|x(/>-(o + [^]}|xoy(/) d9) 

The two summation terms may be viewed as channel measurements. With this time varying 
response, the measured impairment correlation matrix (Step 114) for the current slot is given by 

where SF is the spreading factor of the traffic data. Scaling by 256/SF is used here to get 
absolute impairment levels. Here only the diagonal elements of this matrix are computed. 
[0032] Note that the effect of using interpolated channel estimates as given above is that the 
channel estimate for measuring impairment correlations across a given interval, e.g., a WCDMA 
slot, can change to reflect, for example, changing fading conditions. This method contrasts with 
the alternative approach where impairment correlations (e.g., interference covariance) is 
computed as the difference between each one in a series of samples taken across the slot and 
a channel estimate value that is unchanging for the entire slot. 

[0033] With this embodiment, we may optionally smooth the measured impairment 
correlation matrix before using it to perform parameter fitting for the impairment correlation 
model. Such smoothing can be done with exponential filtering. This can be considered part of 
Step 114. Once this is done, we denote the diagonal element corresponding to finger/as r{f). 
If smoothing is used, then we are basically fitting a smoothed measurement of the impairment 
correlation matrix to a smoothed model of the impairment correlation. Next, we compute the 
fitting parameters. Rather than using LS joint fitting approach, we simply fit (update) the model 
by determining one parameter and then the other. 
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[0034] This method of determining the model fitting parameters exploits the fact that the 
term /? essentially is the noise power, which may be obtained using the diagonal elements of the 

measured impairment correlation matrix, R , according to the method given below. While such 
details provide specific methodology, the broad approach is best understood as obtaining a first 
summation of terms that include both noise and interference powers in first proportional 
amounts and then subtracting a second summation of terms that include noise and interference 
powers but in a different proportion than the first such that subtracting the second summation 
from the first yields a reasonably good estimate for J3. 

[0035] To do this, we partition the set of fingers, denoted F, into two subsets: P and E. The 
set P corresponds to fingers placed on signal paths. The set E corresponds to "extra" fingers 
placed off of the signal path delays. Also, we use \l\ and \e\ to denote the number of fingers in 
these subsets. We compute 



With this embodiment, smoothing to obtain ft is optional — thus processing may use f3 ms t or a 
smoothed version of it. 

[0036] Now, to obtain the second parameter, or, we use 



where the numerator is a function of the difference between the element corresponding to finger 
/in the measured covariance matrix and the corresponding element in the modeled noise 
correlation matrix as weighted by a function of p . The weighting given by q overemphasizes the 
noise by a desired amount and provides an added measure of stability to the problem solution. 
The more weight given, the more the solution leans toward a RAKE solution. Also, the 
denominator is given as w, which is a function that optionally smoothes the argument over time. 




(21) 



H(r{f)-q{p inst )r n (f)) 




(22) 
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Ideally, this should smooth the modeled correlation at the same level as the measured 

impairment correlation is smoothed. 

The mapping function q for example can be a linear function given as 

q(x) = ax + b. (23) 
It is recommended that a be chosen in the range [1 ,2] and b be chosen in the range [0,0.1]. A 
good choice is a = 2, b - 0. Optionally, a insi may be smoothed to obtain a . 
[0037] Another embodiment of the present invention may have particular advantages in the 
context of WCDMA's High Speed Downlink Shared Channels (HS-DSCHs). Thus, Fig. 3 
illustrates exemplary processing logic in Steps 140-156 that largely corresponds to the steps 
given in Fig. 2 with the exception that such operations consider certain state information, e.g., 
Steps 150-154 are based on state-specific information as detailed below. 
[0038] With HS-DSCH, a WCDMA base station allocates a portion of its power to high- 
speed packet data access. For example, it may allocate seventy percent of its power to voice 
users and transmit to them all the time. The remaining thirty percent of its power may be 
allocated to HS-DSCH on an as needed basis. That is, the base station only transmits the HS- 
DSCH signal when there are packets to send. Thus, the base station generally operates in one 
of two states: a full-power state when it is transmitting voice traffic and HS-DSCH traffic, and a 
reduced-power state when it is not transmitting HS-DSCH traffic. 

[0039] The exemplary receiver can be configured to use state-dependent impairment 
correlation model information. A mobile terminal using the HS-DSCH service normally monitors 
control information that tells it which state the base station is in. Thus, the mobile terminal can 
determine which state-based information to use for calculation of the impairment correlation 
model. 

[0040] In an exemplary embodiment, the impairment correlation model can be based on the 
same two impairment terms as given above, i.e., an interference term and a noise term, but the 
model fitting parameter a is different between the two states. Thus, the mobile terminal can 
maintain two separate parameter values, a { and a 2 . When demodulating a given slot or a 
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particular Transmit Time Interval (TTI), the mobile terminal can use the parameter 

corresponding to the state of the base station. 

[0041] Thus, with respect to Fig. 3, when the model fitting parameters are estimated, the 
instantaneous value is used to update (smooth) the corresponding model fitting parameter value 
that corresponds to the state of the system. This updated value is then used to form the model 
H(slot) that is used to generate the RAKE combining weights. Initially, a single value can be 
estimated and then used to initialize the second parameter. Of course, other forms of 
initialization are possible. Similarly, the mobile terminal may wish to keep multiple values of p , 
corresponding to different states of neighboring base stations. Considering state information for 
neighboring base stations has value, for example, where the mobile terminal monitors multiple 
base stations to determine whether it should be served by another base station. 
[0042] For a mobile terminal not using HS-DSCH service, several options are available. 
One is to simply use the first embodiment of the present invention as exemplified by the 
processing logic of Fig. 2. The model fitting parameter values can be averaged over the two 
states of base station operation. Another option is for such mobile terminals to monitor HS- 
DSCH control information, and use the state information to have state-dependent model fitting 
parameters according to the logic of Fig. 3, for example. Yet another option is for the mobile 
terminal blindly to estimate which state the base station is in and maintain multiple model fitting 
parameter values accordingly. The state can be estimated blindly by measuring impairment 
power and thresholding this to form different states, for example. 
[0043] With the above exemplary embodiments in mind, Fig. 4 broadly illustrates a 
transmitter 10, e.g., a wireless network base station, and a receiver 12, e.g., a mobile terminal or 
other wireless communication device. An exemplary transmitter 10 uses spread-spectrum 
modulation to transmit spread-spectrum signals including a traffic signal and a pilot signal. The 
signals pass through a radio channel and are received at one or more antennas of receiver 12. 
Thus, receiver 12 receives the transmitted signal(s) plus noise and interference, and a radio 
processor 14 generates received signal samples, r, from that corrupted received signal. While 

the details of radio processor 14 in at least some embodiments will be familiar to those skilled in 
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the art, an exemplary radio processor 14 comprises filtering and conversion circuits, e.g., 

analog-to-digital converters, such that the received signal is represented by a series of digitized 

baseband signal samples input to a G-RAKE processor 16. In turn, processor 16 demodulates 

the received signal samples to produce soft values or bit estimates. These estimates are 

provided to one or more additional processing circuits 18 for further processing, such as 

forward-error-correction (FEC) decoding and conversion into speech, text, or graphical images, 

etc. Those skilled in the art will recognize that the particular information type(s) carried by the 

received signal and the particular processing steps applied by receiver 12 are a function of its 

intended use and type. Processor 16 also estimates received signal quality by producing SIR 

estimates that also are based on the modeled impairments. 

[0044] Fig. 5 illustrates an exemplary processor 16 comprising a correlation circuit 20, a 
combiner circuit 22, a finger placement circuit 24, and a combining weight and SIR estimate 
generation circuit 26. In operation, received samples are provided to the finger placement circuit 
24, which determines what delays to use in correlating the received signal samples to the traffic 
spreading sequence(s). These delays are used in correlation circuit 20 to delay either the 
received signal sample stream or the spreading codes used by the individual RAKE fingers in 
the correlation circuit 20 to produce traffic correlation values according to the various finger 
placements. The delays are also provided to the combining weight and SIR estimate generation 
circuit 26, which computes the RAKE combining weights used to combine the RAKE finger 
output signals from the correlation circuit 20. It also computes an SIR estimate. The finger 
output signals from correlation circuit 20 are combined using the combining weights in the 
combiner 22, producing combined values or soft bit values. 

[0045] An exemplary correlation circuit 20 comprises a plurality of correlation units, also 
referred to herein as RAKE fingers, and each correlation unit can be placed at a desired relative 
time offset with respect to the received signal using correlation code offsets and/or adjustable 
delay elements, such as buffers. In exemplary G-RAKE operation, finger placement circuit 24 
controls correlation circuit 20 such that one or more RAKE fingers are time aligned with the 
relative path delays of selected signal images in the received signal (on-path fingers) and, 
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typically, one or more of the RAKE fingers are placed off-path. Processor 16 can be configured 
to adjust the alignment of on-path and/or off-path RAKE fingers dynamically to maximize the 
SNR of the RAKE combined signal output from the combiner circuit 22. 
[0046] On that point, G-RAKE processor 16 offers improved performance compared to 
conventional RAKE receivers under at least some reception conditions by considering the 
effects of received signal impairment correlations between the RAKE fingers in its generation of 
the RAKE combining weights. Of course, in the context of the present invention, combining 
weight generation benefits from the use of modeled impairment correlations. To that end, Fig. 6 
illustrates an exemplary combining weight and SIR generation circuit 26 that is configured to 
generate RAKE combining weights according to one or more exemplary embodiments of model- 
based signal impairment compensation. 

[0047] According to the illustration, circuit 26 comprises a correlation circuit 30, a channel 
tracker circuit 32, a SIR calculator 33, a combining weight calculation circuit 34, an impairment 
correlation estimation circuit 36, a structured element calculation circuit 38, a model fitting 
parameter calculation circuit 40, and a modeled impairment calculation circuit 42. These latter 
three elements functionally cooperate as a "impairment modeling circuit" that may be 
implemented in hardware and/or in software. 

[0048] In operation, the received samples are provided to correlation circuit 30, which 
correlates the received samples to a pilot or other reference signal spreading sequence, and 
removes symbol modulation producing pilot correlation values. The channel tracker circuit 32 
receives these pilot correlations and uses them to estimate or otherwise track channel 
coefficients for the received signal. These coefficients are provided to estimation circuit 36, 
which also receives the pilot correlations. Circuit 36 is configured to obtain impairment samples, 
i.e., impairment measurements for the received signal, by subtracting channel estimates from 
the pilot correlations and further configured to calculate impairment correlation measurements 
by correlating the impairment samples with each other and with themselves. 
[0049] Structured element calculation circuit 38 receives the channel estimates and uses 
them to construct the elements corresponding to the impairment terms R 7 and R„ of the 
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impairment model, which are referred to herein as structured elements. The impairment 
correlation measurements, i.e., matrix R , and the structured elements are provided to the model 
fitting parameter calculation circuit 40, which uses them to form the model fitting parameters a 
and p . The fitting parameters and the structured elements are provided to the modeled 

impairment calculation circuit 42, which produces the modeled impairment correlation matrix R . 
The channel estimates and the modeled impairment correlation matrix are provided to the 
weight calculation circuit 34, which produces the combining weights to be used by combiner 
circuit 22 in RAKE combining the finger output signals from correlation circuit 20. The channel 
estimates and the modeled impairment correlation matrix are also provided to the SIR calculator 
33, which produces a SIR estimate for power control purposes. 

[0050] In an exemplary embodiment, a receiver circuit to determine received signal 
impairment correlations for use in received signal processing comprises the impairment 
correlation estimator 36, which in an exemplary configuration measures received signal 
impairment correlations for a received signal of interest, and the one or more impairment 
modeling circuits (e.g., circuits 38, 40, and 42), which in an exemplary configuration implement a 
model of received signal impairment correlations that comprises one or more impairment terms 
scaled by corresponding model fitting parameters. The exemplary impairment modeling circuits 
further adapt each of the model fitting parameters responsive to recurring measurements of the 
received signal impairment correlations as provided by the impairment correlation estimator. As 
illustrated, such a receiver circuit may be included in or associated with a RAKE processor 
and/or with a SIR estimator, both of which may use the modeled received signal impairment 
correlations in their operations. 

[0051] With the above exemplary details in mind, those skilled in the art will appreciate that 

the present invention broadly fits measured received signal impairment correlations, which 

typically have some estimation errors, to an impairment model comprising a structured form that 

is the scaled sum of structured matrix elements. Thus far, the use of two structured matrices 

has been illustrated, one representing an interference term and the other a noise term. In the 

context of a mobile receiver operating in a cellular communication network, the interference term 
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may represent own-cell interference and the noise term may represent white noise and other 
interference. As will be shown below, that model can be extended to include an other-cell 
interference term scaled by a corresponding model fitting parameter. 
[0052] In general, the method can be extended to include additional model terms by 
including other structured elements corresponding to, among other things, other-cell 
interference. In modeling one additional base station, Eq. (4) becomes: 



where R 0 corresponds to other-cell interference. The structure for other-cell interference is 
slightly different than that of own-cell interference, as there is no orthogonal code property for 
this form of interference. As a result, the elements of R 0 are given by, 



where the tilde over the g terms indicates that they are medium channel coefficients 
corresponding to the channel from the other-cell base station to the receiver which has L paths. 
The exemplary receiver can be configured to estimate these terms by correlating to the pilot 
signal of the other-cell base station. Notice also that m=0 is not excluded in the last summation. 
[0053] In Eq. (25) above, it is assumed that the medium channel coefficients corresponding 
to the channel from the other-cell base station to the receiver are estimated. In cases where the 
instantaneous medium channel coefficients are not available, R 0 (d l9 d 2 ) can be formed through 
any number of alternatives, including the following exemplary alternatives. 
[0054] If the average path strength and delays corresponding to the channel from the other- 
cell base station to the receiver are available from the searcher, Ro(d u d 2 ) can be formulated as 



[0055] If the path searcher only provides the delay corresponding to the dominant echo, 
R Q {d u d 2 ) can be formulated as 



R « oRj + j3R n + jR 0 



(24) 



*o (di ,d 2 )= II?/?/ I *p id, - mT c -?,)*/ (*2 -mT c - r J , 



(25) 




(25a) 



RoWd 2 )= f R p {d x -mT c -T Q )R/(d 2 -mT c -? 0 ). 



(25b) 
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In this case, the path strength is folded into 7 . 

[0056] Another way to extend the method to other-cell interference is to model the other-cell 
interference as white noise that has passed through the transmit pulse shaping filter. This 
approach does not require estimating the channel response of another base station. The 
resulting R 0 (d Xi d 2 ) is simply: 

R 0 (d u d 2 ) = R q (d x -d 2 ) f (26) 

where R q (x) is the autocorrelation function of the pulse shape convolved with itself, possibly 

normalized so that the zero-lag element is one. 

[0057] Note that it is possible to utilize a combination of these different approaches. In this 
case, the correlation matrix of other-cell interference can be expressed as a weighted sum of the 
^0(^1^2) in equations (25), (25a), (25b) and (26). 

[0058] When the exemplary receiver is in soft handoff between two or more base stations 
there are two or more received signals of interest, e.g., traffic signals, being transmitted for the 
receiver by the network from different radio sectors or from different sites. 
[0059] In such cases the receiver assigns one set of its RAKE fingers to extract one of the 
signals and another set of its RAKE fingers to extract the other. These sets of fingers can be 
treated separately, as far as combining weight computation is concerned. Thus, the present 
invention can be applied separately to each finger set. 

[0060] In such soft handoff cases, the receiver performs channel estimation for multiple 
base stations. Thus, the receiver could be configured to include other-cell interference in its 
impairment model, e.g., include the effects of other-cell interference in the model's impairment 
correlation matrix. When computing combining weights for the first base station signal, the 
receiver would treat the second base station signal as other-cell interference. When computing 
weights for the second base station signal, it would treat the first base station signal as other-cell 
interference. 

[0061] Besides soft handoff, another way that the receiver receives multiple transmitted 
signals is when transmit diversity is used. In WCDMA, there are basically two forms of transmit 
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diversity. One form is open loop or STTD, which uses an Alamouti approach to code 
information across the two transmit antennas. Another form is closed loop, which relies on 
feedback so that the same symbols sent from the two antennas arrive in-phase at the receiver. 
In either case, there are effectively only 5 pilot symbols instead of 10 in a slot for each transmit 

antenna. Thus, in equation 1 1 the summation would be from 0 to 4 and would be replaced 
by -. For Eq. 19, we can use 



[0062] Also, in either case, the impairment correlation matrix is modeled like in the multiple 
base station case, with a term for each transmit antenna (see (24)). However, in (25), the m-0 
term should be excluded, just as in (15). For additional exemplary details on these and related 
calculations, one may refer to the co-pending and commonly assigned U.S. patent application, 
which is titled "A METHOD AND APPARATUS FOR RECEIVED SIGNAL QUALITY 
ESTIMATION," and which is assigned Attorney Docket No. 4015-5191 . That application, filed 
on even date herewith, is incorporated herein by reference in its entirety. 
[0063] One approach is to have only one R matrix of the form of (24). Details of how this 
matrix is used are given below. Note that this model can be used whether or not the particular 
receiver is using transmit diversity or not. Thus, as long as the base station is transmitting 
something on two antennas, then such an extended model of the impairment correlation matrix 
should be used. For the STTD approach, data symbols are transmitted in pairs. On transmit 
antenna 1 , symbol 1 is sent in the first symbol period, whereas symbol 2 is sent in the second 
symbol period. On transmit antenna 2, the negative conjugate of symbol 2 is sent in the first 

symbol period, whereas the conjugate of symbol 1 is sent in the second symbol period. Let 

be the response from transmit antenna A and h B be the response from transmit antenna B. 
Then the following two weight vectors are formed: 

w A =R- l (s!ot)h A (28) 



5 




(27) 
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w fl = K~ l (slot)h B . (29) 

Let y, and y 2 be the vectors of traffic despread values during the first and second symbol 

periods. Then the decision variables for symbols 1 and 2 are obtained by 

*i=w? yi +(w?y 2 )T (30) 

*2=wZy 2 -(w2yJ. (31) 
Thus, the impairment correlation matrix would be used to form two combining weight vectors. 
The SIR can be estimated using 

SlR = h A H ii- x {slot)h A + h^R" 1 {slot% B =h A N w A + h B H w B . (32) 
Thus, the SIR would be the sum of two SIR terms. For the closed loop approach, the same 
symbol is transmitted using the same spreading code from the two transmit antennas. In one 
closed loop approach, the relative phase of the two transmissions is adapted. In the other 
approach, relative amplitude is also changed. For simplicity, we can think of this as sending the 
symbol s on one antenna and symbol es on the other where e is a complex quantity. The 
receiver would know e or could learn it from the transmissions. From the two pilot channels, the 
receiver would estimate two responses and h B . The combining weights would then be 
formed using 

w = R" 1 {slot ifc A + eh B )= R" 1 {slot)h x . (33) 
SIR would be estimated using 

SIR = h^R" 1 (slot)h x = h x H w . (34) 
[0064] In general, with transmit diversity we have a number of options for computing the 
impairment correlations. Broadly, the exemplary receiver may be configured in one or more 
ways depending on the diversity scenario. Where there are separate (RAKE) finger locations for 
each transmit diversity signal, we may use separate impairment correlation models with 
separate fitting. Conversely, where the finger locations are the same for all transmit diversity 
signals, we can "stack" Eq. (6) by adding elements to the vectors and matrices accounting for all 

transmit diversity signals. Pre-combining the impairment correlation measurements, 
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Ri ,R 2 , ,R* i for transmit diversity signals 1 ..n represents another approach. For the in- 
between case wherein only some fingers positions between the transmit diversity signals are the 
same, the exemplary receiver can be configured to use subsets of a correlation matrix formed 
based on the fingers used for each transmit diversity signal. 

[0065] For the separate modeling case, the exemplary receiver models the impairment 
correlations separately for each transmit antenna. Thus, using transmit antennas 1 and 2 as an 
example, the receiver would determine the impairment correlations as 

R x =a\R/ x +fi\R n + Xi R /2 for Antenna 1, and as R 2 = <*2 R /i +/ ? 2 R /i + /2 R /2 for Antenna 2. 
Using per-antenna models may be particularly advantageous where the transmit antennas are 
physically separated from each other. Making the necessary impairment correlation 
measurements needed for model fitting may be based on receiving separate pilot signals from 
the different antennas. 

[0066] For cases where there is only a partial overlap of finger positions between the 
different antennas, the receiver may form the combining weight vectors based on subsets of an 
impairment correlation matrix. For example, assuming that receiver finger positions 0, 1 , and 2 
were assigned to Transmit Antenna 1 and finger positions 0, 1 , and 3 were assigned to Transmit 
Antenna 2, we see that positions 0 and 1 overlap between the two antennas but positions 2 and 
3 do not. In this instance, the impairment correlation matrix R would include four columns (0, 1, 
2, 3) and four rows (0, 1,2,3) but only the corresponding 3x3 subset of that matrix would be 
used for each antenna to generate the combining weights and SIR estimates. 
[0067] If the receiver includes multiple receive antennas, it can be configured to assign one 
or more of its RAKE fingers to one receive antenna signal, and assign one or more remaining 
fingers to each of the other receive antennas. In other words, the exemplary receiver can 
allocate a subset of its RAKE fingers to each receive antenna, and can be configured to perform 
channel estimation separately for each antenna. In such cases, impairment correlation 
measurements can proceed as before except that some of the correlation measurements 
correspond to correlations of impairment between different receive antennas. 
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[0068] Similarly, the overall channel estimates comprise sets of channel estimates 
corresponding to the different receive antennas. Thus, let h t and h 2 represent the net channel 



impairment correlation matrix between fingers on antenna i and j branches, i.e., the (m,n)-\h 
element of R. . is the interference correlation of RAKE finger m on the j antenna branch and 

finger n on the j antenna branch. In this case, the RAKE combining weights can still be 
formulated as 



[0069] Impairment correlations among fingers from the same antenna branch ( R n and 
R 22 ) can be estimated using methods discussed thus far. Impairment correlation between 
fingers from different antenna branches can be estimated as follows. 
[0070] First a measurement of R 12 is computed from impairment realizations, which 

measurement may be denoted as R 12 - From the foregoing information, it should be understood 
that R l2 can be expressed as a weighted sum of an own-cell interference component, R I2 / and 
an other-cell interference component, Ri2,o> 



[0071] The weighting factor oris determined by own-cell power, whereas 7 is determined by 
other-cell interference power and, in some cases, also by the correlation between antennas. 
The white noise can be assumed to be uncorrelated between antennas. Note that the same a 
and y appear in all four sub-matrices. When antenna gains differ, it may help to have different 
p values for R n and R 22 . 

[0072] The (ij)-th element of R 12>/ and K\ 2i o can be computed from net response and 
pulse autocorrelation as 



responses for first and second receive antennas, respectively. Further, letR, , be the 



w = R" l h, 



(35) 




R 12 -oR 12>/ + jRi2,o- 



(36) 
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R\20{di>dj)= Z Zgijglg llRpidi-mT, -mT c -r q ) 



(37) 



where g it is the medium channel coefficients corresponding to the /lh path from the own-cell 

base station to the receiver antenna /, and g it , is the medium channel coefficients 

corresponding to the /th path from the other-cell base station to the receiver antenna L 
[0073] Note that in Eq. (37), it is assumed that the medium channel coefficients 
corresponding to the channel from the other-cell base station to the receiver are estimated. In 
cases where the instantaneous medium channel coefficients are not available, R l2i o[ d i> d j) can 

be formed by other means, such as through one of the following alternatives. 
[0074] If the path average strength and delays corresponding to the channel from the other- 
cell base station to the receiver are available from a searcher circuit in the receiver, then 
R\isX d i^ d j) can be formulated as 



In this case, the correlation between antennas is folded into 7 . 

[0075] If the path searcher only provides the delay corresponding to the dominant echo, 
R\2,o[ d i> d j) can be formulated as 



In this case, both the correlation between antennas and the path strength are folded into 7 . 
[0076] Another way to extend the impairment model to include impairment correlations 
arising from other-cell interference is to model the other-cell interference as white noise that has 
passed through the transmit pulse shaping filter. This approach does not require estimating the 




(38) 



*\wW*j)= HRpfa -mT c -T 0 ]R p *{dj -mT c -f 0 ) . 



(39) 



channel response of another base station. The resulting R^o^i^j) is simply: 
RnA d i> d j)= R M- d j)* 



(40) 



26 



Ericsson Ref. No. P18204-US1/P 19093 
Coats & Bennett Docket No. 4015-5161 

[0077] Note that it is possible to utilize a combination of these different approaches. In such 
cases, the correlation matrix of other-cell interference can be expressed as a weighted sum of 
the R\2 t o[ d i> d j ) in E Q S - (37), (38), (39) and (40). Further, if finger delays are spaced no closer 
than three-fourths of a chip apart, the other cell interference covariance can be approximated as 
R i2,o = 1 • WithR 12 , R 12/ and Ri2,a » unknown weighting factors orand j^can be solved using 

the least squares approach. Preferably, R n , R 12 , and R 22 are used together to solve for a , 
J3 , and 7 . 

[0078] Of course, these multi-transmitter and multi-antenna embodiments represent just 
some of the many possible variations that can be practiced in accordance with the present 
invention. Those skilled in the art will recognize additional features and advantages falling 
within the scope of the present invention. Fundamentally, the present invention comprises a 
method and apparatus wherein received signal impairment correlations are compensated for in 
RAKE combining weight generation using a model-based approach that yields computational 
and performance advantages. Such advantages are gained at least in part by forming the 
model as one or more structured elements scaled by model fitting parameters allowing the 
model to be efficiently adapted at each of one or more successive time instants based on fitting 
the parameters in accordance with measured impairments. 

[0079] The impairment model may be configured to consider multiple sources of 
interference, such as same-cell and other-cell interference, and may be configured for operation 
with respect to multiple transmitted signals from two or more radio sectors and/or base stations 
in a cellular radio environment, and/or with respect to multiple transmitted signals as obtained 
from two or more receive antennas. Fig. 7 provides an exemplary illustration of a wireless 
communication network 50, which may be configured as a WCDMA wireless cellular network, as 
an IS-95/IS-2000 wireless cellular network, or configured according to some other open or 
proprietary communication standard. 

[0080] Network 50 supports a plurality of mobile terminals 52 — only one terminal is shown 
for clarity — operating in one or more radio service areas denoted as C1 , S1 , C2, S2, and so on, 
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to represent Cell 1 , Sector 1 and Cell 2, Sector 2, etc. Those skilled in the art will appreciate 

that the terms "cell" and/or "sector" as used herein should be given broad construction and, in 

general, the term sector should be understood as identifying a given radio coverage area at a 

given carrier frequency. Thus, a given cell may have multiple overlaid radio sectors 

corresponding to multiple carrier frequencies. 

[0081] In any case, network 50 communicatively couples mobile terminals 52 to one or more 
external networks 54, such as the Public Switched Telephone Network (PSTN), the Internet or 
other Public Data Networks, ISDN-based networks, etc. Such coupling is supported by Radio 
Access Network (RAN) 56, which provides the radio link(s) to the mobile terminals 52, and 
which interfaces to one or more Core Networks (CNs) 58 that in turn link to the external 
networks 54. Those skilled in the art will appreciate that the particular network architecture 
adopted, and the entity nomenclature used varies according to the network standard involved, 
but such variations are not germane to understanding or explaining the present invention. 
Moreover, it should be understood that the illustrated network is simplified and that actual 
network implementations likely will have additional entities not illustrated herein for clarity. 
[0082] Nonetheless, an exemplary RAN 56 comprises one or more base station systems, 
each typically comprising a control entity and one or more distributed radio transceiver entities. 
In the illustration, such entities are depicted as one Base Station Controller (BSC) 60, and a 
plurality of associated Radio Base Stations (RBSs) 62, e.g., 62-1 , 62-2, and 62-3. The 
illustrated mobile terminal 52 includes the exemplary receiver 12 illustrated earlier herein, which 
it may implement using various processing circuits, including A/D converters, filters, DSPs or 
other digital processors, memory, and the like. In at least one exemplary embodiment, mobile 
terminal 52 includes one or more DSPs and/or Application Specific Integrated Circuits (ASICs), 
or other programmable devices, to implement receiver 12, including G-RAKE processor 16, as 
illustrated in exemplary fashion in Figs. 4, 5, and 6. It should be understood that at least a 
portion of the present invention's functionality thus can be embodied as stored computer 
instructions in the form of micro-code, firmware, software, etc. 
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[0083] More generally, the present invention can be implemented in hardware, software, or 

essentially any combination thereof, according to the needs of particular design. Indeed, the 

present invention is not limited by the foregoing discussion or by the accompanying figures. 

Instead, the present invention is limited only by the following claims and the reasonable 

equivalents thereof. 
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